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This paper defines a notion of parallel transport in a lattice of quantum par- 
ticles, such that the transformation associated with each link of the lattice is 
determined by the quantum state of the two particles joined by that link. We 
CS| ■ focus particularly on a one-dimensional lattice — a ring — of entangled rebits, 

which are binary quantum objects confined to a real state space. We consider 
states of the ring that maximize the correlation between nearest neighbors, 
and show that some correlation must be sacrificed in order to have non-trivial 
parallel transport around the ring. An analogy is made with lattice gauge 
theory, in which non-trivial parallel transport around closed loops is associ- 
q-i ated with a reduction in the probability of the field configuration. We discuss 

the possibility of extending our result to qubits and to higher dimensional 
lattices. 
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1 Introduction 



In lattice gauge theory, the gauge field assigns to every pair of neighboring 
lattice sites a transformation that tells how to "transport" a vector, repre- 
senting an internal property such as quark color, from one site to the other. 
That is, with every ordered pair (j, k) of neighboring sites, or link, one as- 
sociates a transformation U(k,j) which is an element of the gauge group. 
Indeed, this assignment of transformations to links constitutes the configura- 
tion of the gauge field. If we apply U(k,j) to a vector v associated with site 
j, we can interpret the image vector, v' = U(k,j)v, as the result of moving 
v from site j to site /cfj This process is called parallel transport, and the 
transformation U(k,j) is sometimes called a parallel transporter. 

This paper is not about lattice gauge theory but about actual lattices 
consisting of simple quantum particles, in which each particle is correlated 
with its nearest neighbors. As we will see below, for typical states of the 
lattice, one can use the state itself to specify a notion of parallel transport. 
Our conception is distinct from the more familiar notion of quantum parallel 
transport (expressible in terms of gauge fields JTJ, [§), in which a quantum 
particle is physically moved either in actual space or in a parameter space 
[0. p3, ||. In our approach there is no physical motion or evolution; rather, 
the transformation that we associate with a link expresses something about 
the relationship between the particles joined by that link. The present work 
is an initial exploration into a possible analogy between the quantum state 
of a lattice of quantum particles and the configuration of the gauge field in a 
lattice gauge theory. Roughly, the analogy we are looking for would be along 
the following lines (it will be spelled out more precisely in later sections of 
the paper). 

A lattice gauge theory assigns a probability distribution to the set of pos- 
sible field configurations, the probability density of a configuration A being 
proportional to exp[— 5(^4)] where the real function S(A) is the action. The 
action depends on the effect of parallel transport around each of the elemen- 
tary plaquettes of the lattice; that is, it depends on the transformations that 
the field configuration associates with these elementary closed paths. The 
more these transformations differ from the identity, the higher the action, 

It helps to read the arguments of U from right to left. This notation makes sense when 
combined with the usual notation for a sequence of operations, in which the operator on 
the right acts first. 
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and therefore the lower the probability of that particular field configuration. 
Let us call this effect "the probability cost of twisting." I am looking for 
something similar in a lattice of correlated quantum objects, but instead of 
a probability cost, I am looking for a "correlation cost." We will be focusing 
our attention on states in which nearest neighbors are maximally correlated. 
The question is whether this maximal correlation must be reduced if there 
is to be non-trivial parallel transport around closed loops (that is, transport 
whose net effect is not null), and if so, whether the reduction in correlation 
becomes more severe as the effect of such parallel transport differs further 
from the identity. In other words, in an entangled lattice, is there a correla- 
tion cost of twisting? 

In fact this paper only begins to answer this question. Although we set 
up the problem for lattices of arbitrary dimension, the one concrete example 
we work out in detail is the case of a one- dimensional ring. Moreover, for 
most of the paper we will restrict our attention to the simplest possible 
quantum object, namely, a rebit, a fictitious object whose state space is a 
two-dimensional vector space over the real numbers 0. Thus we will mostly 
be analyzing a closed ring of rebits. At the end of the paper I discuss the 
generalization to qubits and to higher-dimensional lattices. 

This research is related to a recent line of work on "entanglement shar- 
ing," concerning the ways in which quantum entanglement can be shared 
among several objects. A number of authors have found constraints on the 
sharing of entanglement that follow directly from the structure of quantum 
mechanics itself and not from any particular Hamiltonian. For example, it 
has been shown that any entanglement that might exist between a pair of 
qubits limits the extent to which either of them can be entangled with a 
third qubit [0, §]. There are similar limits for n qubits all entangled with 
each other || [10], [11]]. Another example is an entangled ring: in a translation- 



ally invariant state of a ring of qubits, there is a certain maximum possible 
entanglement between nearest neighbors [JT2]. (The measure of entanglement 
used in all of these studies is the entanglement of formation 0].) In 
this paper we are putting a different sort of condition on the quantum state 
of a multipartite system — we are imposing a certain configuration of twists 
in the nearest-neighbor correlations — and we are asking what constraint this 
condition places on the strength of the correlations. 

I hope that the results of this research will ultimately be useful in analyz- 
ing systems of entangled particles on a lattice, such as magnetic systems. If 
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there are simple laws of quantum correlations that transcend any particular 
Hamiltonian, then these laws might lead to the identification of interesting 
generic properties of quantum many-body systems. Arguments along these 
lines, particularly focusing on entanglement, have appeared recently in the 



literature |L5] , jig , [17], |i8fl . But at least as much of the actual motivation for 
the present work comes from pure curiosity: I wonder how close an anal- 
ogy one can draw between the degree of correlation between particles in a 
quantum lattice and the probability density of a field configuration in lattice 
gauge theory. Of course there are many connections between lattice gauge 
theory and the theory of many-body systems — see, for example, Ref. |]19[ and 
references cited therein — but I am looking for an analogy along the particular 
lines traced out above. 

The reader may have noticed that in describing the work to be presented 
here (as opposed to earlier work), I have been using the word "correlation" 
rather than "entanglement." Though they are related, the two concepts 
are not the same. In this paper I focus on correlation because it seems 
natural in this context and it is easy to work with. But it would also be 
interesting to explore the same questions using one of the standard measures 
of entanglement. I might add that the states we will primarily be concerned 
with are in fact highly entangled; hence the reference to an "entangled ring" 
in the title. 

The paper is organized as follows. First we review briefly those aspects 
of lattice gauge theory that have suggested our main question. We then 
define a rebit more precisely and develop our notion of parallel transport. 
We analyze in some detail the case of a ring of rebits and determine whether 
non-trivial parallel transport does indeed entail a "correlation cost." Finally 
we ask how the problem and the results are likely to change when extended 
to more complex systems. 



2 Lattice gauge theory and a simple analogy 

Ideally, lattice gauge theory is done on a four-dimensional lattice represent- 
ing spacetime, except that the fourth dimension represents imaginary time, 
so that it acts in many respects like another spatial dimension. The results 
of a calculation can be interpreted in terms of real time by means of analytic 
continuation. One consequence of the use of imaginary time is this: in com- 
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puting the expectation value of an observable, one does not sum up complex 
amplitudes associated with different histories; rather, one takes a weighted 
average of the observable of interest using real weights ||. As mentioned 
above, the weighting function is proportional to exp(— S), where the action 
S is a real function of the field configuration. More precisely, in a pure gauge 
theory, in which there are no matter fields but only the gauge field itself, the 
expectation value of an observable B is 

(B) = ^fBe- 8 ^UdU(k 7 j). (1) 

<J,k) 

Here A is the configuration of the gauge field, which assigns a parallel trans- 
porter U(k,j) to each link (j, k). Each such transformation U is an element 
of the gauge group, e.g., U(l) for electrodynamics or SU(3) for chromody- 
namics. The integral in Eq. ([TJ) is over all field configurations, that is, over all 
possible parallel transporters for each link, and dU(k,j) indicates the invari- 
ant measure over the gauge group. In the integral, the parallel transporters 
for different links are independent, except that U(k,j) = U(j, A;) -1 , so that 
only one of these two ordered pairs needs to be represented in the integral. 
The normalizing constant Z is simply 

Z= fe- s ^l[dU(k,j). (2) 

The action S depends only on the results of parallel transport around 
plaquettes of the lattice, e.g., elementary squares in a cubic lattice. Let V be 
a plaquette, which we can think of as a sequence of lattice sites; for a cubic 
lattice V would consist of four sites j, k,l,m. Let U(V) be the net effect of 
parallel transport around plaquette V; in the square example, U(V) would 
be U(j,m)U(m,l)U(l,k)U(k,j) (the right-most operator acting first). S is 
defined so that it increases as the plaquette transformations U (V) get farther 
from the identity. Various action functions with this property have been used 
in the literature; the one originally proposed by Wilson for an SU(iV) gauge 
field is @] 

Soc-£Re(Tr([/(7>))), (3) 
v 

where the sum is over all plaquettes in the lattice. Note that if V consists 
of the same set of points as V , but with a different starting point or with 
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the points taken in the opposite order, then Ke(Tr(U(V'))) = Re(Tr(£/(P))). 
Thus we need include in the above sum only one ordered set V representing 
each geometric plaquette. 

A gauge transformation associates with each lattice point j a group ele- 
ment A(j), and under such a transformation each parallel transporter U(k,j) 
transforms according to 

U(k,j)^A(k)U(k,j)A(j)-\ (4) 

It is easy to see that, though a gauge transformation changes the field con- 
figuration, it does not change any of the plaquette transformations U(V). 
Therefore it does not affect the action S and so has no physical consequences. 

This invariance under gauge transformations provides a simple analogy 
between the configuration of a gauge field and the state of a lattice of quan- 
tum particles. Consider, for example, a lattice of qubits. Rotating each of 
the individual qubits separately is analogous to a gauge transformation. The 
state of the lattice changes under such rotations, but certain physical prop- 
erties do not change. In particular, any reasonable measure of the degree of 
entanglement or correlation between two qubits does not change. So at least 
in this one modest respect, the degree of correlation in a quantum lattice is 
similar to the action or the probability density of a field configuration in a 
lattice gauge theory. We want to see whether the similarity goes any further 
than this. 

3 What is a rebit? 

As we have said, the quantum object we will mostly be concerned with in 
this paper is the rebit. We now define this object more precisely. 

A pure state \ip) of a single rebit is simply a normalized vector in a two- 
dimensional real vector space. A mixed state of a rebit is a mixture of pure 
states: 

P = ^Pi\A)(A\, (5) 

i 

where Pi > and J2iPi — 1- Equivalently, a mixed state can be represented 
as a real, symmetric 2x2 matrix with unit trace and no negative eigenvalues. 

Of course any rebit state is also a qubit state, and for our purposes it will 
be helpful to think of rebits simply as restricted qubits. On the Bloch sphere, 
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the restriction to real density matrices becomes a restriction to the x-z plane. 
But it will be more useful to change the representation by rotating the Bloch 
sphere. Let us rotate all states by 90° in the left-handed sense around the 
positive x axis, so that our rebit states now lie in the x-y plane. A general 
mixed state lying in this plane can be written as 

p= l[l + aa x + ba y ], (6) 

where the cr's are Pauli matrices and the real numbers a and b satisfy 
a 2 + b 2 < 1. Let us call this representation of rebit states the "horizon- 
tal representation," as opposed to the original "real-number representation." 
Given a rebit density matrix p expressed in the horizontal representation, we 
can always re-express it as a real density matrix p rea i simply by reversing the 
rotation around the x axis: 

Preal = U pU\ (7) 

where 

Here we have written the matrix in the basis of eigenstates of a z , {|j), |!)}. 
Note that Ua y W = a z and Ua x W = a x . 

In the real-number representation, a pure state of n rebits is a real vector 
in 2 n dimensions, and a mixed state of such a system is a real, symmetric 
2 n x 2 n density matrix. But again we will usually work in the horizontal 
representation, in which each rebit has been rotated by 90° around the x axis. 
It will be helpful to have a simple way of recognizing whether a given n-qubit 
state is a legitimate n-rebit state expressed in the horizontal representation. 
Conceptually, the test is straightforward: apply the above transformation 
U to each qubit — that is, rotate each qubit by 90° (in the right-handed 
sense) around the positive x axis — and see whether the resulting state has 
only real components. Thus for a pure state I*), we insist that [/® n |\&) be 
real in the standard up-down basis. But this is the same as saying that 
U® n \ty) = (£7 _1 )® ri |\I'*) , where the asterisk indicates complex conjugation 
in the standard basis. Multiplying both sides of this equation by U® n and 
noting that U 2 = —icr x , we arrive at the following criterion: 

(-1(7,)®"]*) = I**)- (9) 
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In practice it will be simplest if we also allow ourselves to use state vectors 
of the form exp(ia)|\l/), where a is real and |\?) satisfies Eq. (H). Though 
such state vectors do not become real when they are transformed by U® n , 
their density matrices do become real. Allowing this possibility leads to the 
following weaker condition on an n-qubit state |\&). 



j3 being any real phase. Note that the matrix a x simply interchanges ||) and 
||). Thus we can recognize a pure n-qubit state |\&) as a legitimate horizontal 
representation of an n-rebit state by checking to see that the coefficient of 
each basis state, e.g., |Ttlt), is the complex conjugate of the coefficient of 
the opposite state, in this case ||J,tl); multiplied by a phase factor that is the 
same for all basis states. Let us call Eq. (|l0|) the "rebit condition" for pure 
states. The corresponding test for mixed states can be obtained by a similar 
argument; one finds that a density matrix p of n qubits is the horizontal 
representation of a legitimate rz-rebit state if and only if 



the complex conjugation again being in the standard basis. 

We conclude this section with a word about rotations of a rebit. View- 
ing the states of a rebit as qubit states confined to the equatorial plane 
of the Bloch sphere, we could take as the allowed rotations all the unitary 
transformations that represent rotations around the z axis, that is, all trans- 
formations of the form 



However, our definition of parallel transport will not be able to distinguish 
unitary transformations that are different only by an overall phase factor; 
so we will call such transformations identical. For definiteness we pick a 
standard representative from each of the resulting equivalence classes: a 
rotation by an angle £ around the z axis, with < £ < 2tt, will be represented 
by the matrix 



af n \q) = e i/3 |#*), 



(10) 




.(gin _ * 
x r 3 



(11) 




(12) 




(13) 
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4 Parallel transport in a lattice of rebits 



We now consider a lattice of rebits, on which we want to define a notion 
of parallel transport. For now the structure of the lattice does not matter, 
as long as any two rebits are identified as either being neighbors or not. 
Thus the lattice is simply a graph. The mathematical object that is to be 
parallel transported is a pure state of a single rebit, which we can picture as 
a direction in the horizontal plane. The parallel transporter associated with 
a link in the lattice will be a rebit rotation — represented in the form flTB|)— 
which we can identify with an element of U(l). We want the assignment of 
parallel transporters to links to be determined by the quantum state of the 
lattice itself. 

We begin with the following scenario. Consider a pair of neighboring 
rebits labeled j and k. They are in some joint state p^ k \ a 4 x 4 density ma- 
trix obtained from the state of the whole lattice by tracing over all the other 
particles. Imagine performing an arbitrary orthogonal measurement on parti- 
cle j. In the standard von Neumann model this measurement brings particle 
j into a certain pure state — one of the eigenstates of the measurement — and 
it also brings particle k into some state, typically a mixed state, determined 
as follows. Let be the eigenstate into which particle j is brought by 

the measurement. Then the post-measurement state of particle k is the 2x2 
density matrix 

w ( fc ) = ( 14 ) 

where P = ( , ip^\TrkP^ k ^\il>^} is the probability with which that particular 
outcome occurs. In Eq. fll4|) the matrix operations are done only in the 
space of particle j. To be more explicit, if r = 0, 1 and s = 0, 1 are indices 
associated with particles j and k respectively, we can write the components 
of lu^ in terms of the components of \ip^) and p^ k ^ as 

^ = ^EvFp!Sv#- (15) 

r,r' 

We want to use this measurement scenario to associate with the two- 
particle state simple rotation U(k,j). First, let M. (for "measure- 
ment") be the mapping defined by Eq. ([14]) , which takes each pure state of 
particle j for which P ^ into a pure or mixed state of particle k; that is, 
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M(ip ij) ) = u; (fe) . If P = 0, let us say for definiteness that M{ip {j) ) = Trjp {jk \ 
though it will not actually matter in what follows. For an arbitrary rebit 
rotation R, we define a function F(R) by 



F(R) 



J(ip\R^M(jj)R\ip)PdiP 



(16) 



Here P is the probability given above, and dtp indicates the uniform measure 
over the circle of pure rebit states, normalized so that J dip = 1. Thus F(R) 
is an average fidelity of A4(ip), not with respect to \ip) itself but with respect 
to a rotated version of \ip). As we will see in the next paragraph, depending 
on the density matrix p^ k \ one of the following two conditions will hold: 
(i) F(R) is independent of R, or (ii) there is a unique rotation R = U that 
maximizes F(R). In case (i), we say that there is no correlation between 
particles j and k. In case (ii), we take the special rotation U that maximizes 
F(R) to be the parallel transporter associated with the link (j,k). In a 
certain sense, U is the rotation that most closely approximates the action of 
M. 



Combining Eqs. fll4j) and (|Tg) and the definition of P, we have 



F(R) 



/(V>|Tr fcP 0-*)|V>)# 



(17) 



It is not hard to show that the denominator is always 1/2. Also, in the 
horizontal representation, the integral in the numerator involves rotating \ip) 
around the z axis, so that we can rewrite Eq. ( |TTD as 



F(R) = 2((^o| ® (V'ol^) 



2vr Jo 



2- 



(M,>®i#o», ( 18 ) 



where 5, = (<#) + of))/2 and |^o) 

is some fixed reference state which for 
definiteness we take to be \ipo) = (|T)+ Now, a general two-particle 

density matrix satisfying the rebit condition (p"T|) is of the form 



.r? 



^2 ^3 \ 
C X 2 
C* 6 Xi 

a / 



6 



..if. n 



(19) 



.7; 
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the representation being in the standard basis {|Tt)j |tl)> |IT)> III)}- The 
average over 7 in Eq. flT8|) has the effect of replacing the Xj's in p^ k ' with 
zero and leaving the matrix elements a, b, and c unchanged. To evaluate 
F(R), we write R explicitly as a rotation around the z axis by some angle £: 

*=(J£)- < 2 °> 

Inserting this matrix into Eq. ([18]) we find that 

F(i2) = i[l + 2|c|cos(£-0)], (21) 

where is the phase of the matrix element c; that is, c = |c|e*^. If c is zero, 
then we find ourselves in case (i) mentioned above: F(R) is independent of 
R. Otherwise F is maximized when the angle of rotation £ is equal to the 
phase <p. According to our prescription, then, the parallel transporter U(k,j) 
is the rotation given in Eq. (|20| ) with £ = <fi. 

We will take as our measure of the degree of correlation between particles 
j and k the quantity 2|c|, which ranges from to 1. Eq. ([H]) makes it clear 
that 2\c\ measures the degree of angular correlation between the two particles. 
We can also interpret this quantity in terms of the more standard correlation 
matrix t^ u = Tr[p^ k \a Mi <8> ov)], in which both p and v take as values the 
axis labels x and y. One finds that (2|c|) 2 = [{t xx + r yy ) 2 + {r xy — T yx ) 2 ]/A = 
[2detr + Tr(r T r)]/4. 

Note that the single complex number c determines both the parallel trans- 
porter (through its phase) and the degree of correlation (through its mag- 
nitude). Typically we will be trying to maximize the magnitude of c for a 
given value of its phase. The fact that it is possible for a link (j, k) to have 
an undefined parallel transporter U(k,j) (case (i) above, where c = 0) will 
not cause any difficulties for the problem we will be studying. We will be 
considering the set of all states that are consistent with a given specification 
of the parallel transporters, i.e., the phases of the c's. If a link has c = 0, 
then we simply say that that link is consistent with any specified phase. 

Our notion of parallel transport has a particularly simple interpretation 
if pV k ' has the form ( |19D with all the x^s equal to zero. This will happen, 
for example, if the state of the lattice is invariant under identical rotations 
of all the rebits. If p^ k ^ has this form, then Eq. ([TJ) yields 

= P U(k,j)\^)(^\U(k,jy + (l-p)(I/2), (22) 
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where / is the 2x2 identity and p = 2|c|. Thus the post-measurement state 
of particle k is simply a rotated and partially depolarized version of 
and the weight p of the pure rotated state is our measure of correlation. 

Let us now imagine transporting, mathematically, a rebit state around a 
closed loop in accordance with the above prescription. As we will see, the fi- 
nal state in such a process need not be the same as the initial state but could 
be rotated by some angle 0. One might ask: How does one interpret physi- 
cally this process of transport, and what is the meaning of the rotation angle 
91 I regard our concept of parallel transport primarily as a mathematical 
notion; nothing is being physically transported. However, in the special case 
considered above, in which the state of the lattice is rotationally invariant, 
one can extract from our definition a simple physical interpretation of the 
net rotation angle. Consider a closed loop of n lattice sites jo,ji, ■ ■ ■ ,j n -i- 
At each of the sites jo,---,j n -2, that is, at all but the last site, perform 
an orthogonal measurement, with outcomes labeled '0' and '1', choosing the 
measurement at j m so that it maximizes the probability of getting the same 
outcome (0 or 1) as at site j m -i- That is, we are trying to minimize the 
expected number of flips from to 1 or from 1 to as we go around the 
loop. Now, at the last site, j n -i, one is faced with a dilemma: there will be a 
measurement that maximizes the probability of agreement between j„-2 and 
jn-i, and there will be a (possibly different) measurement that maximizes 
the probability of agreement between j n _i and jo- The angle between these 
two measurements, that is, between their '0' eigenstates, is the angle 8 asso- 
ciated with parallel transport around the loop. We can think of this angle 
as measuring the net "twist" in the nearest-neighbor correlations. 

When the state of the lattice is not rotationally invariant, the angle be- 
tween the two competing optimal measurements at the last site (the two 
being reckoned optimal from different directions) may depend on the choice 
of the initial measurement on the first particle. So the interpretation in this 
case is not as simple. Still, it is reasonable to think of the net rotation an- 
gle as a measure of the net twist in the correlations. In the following section 
where we analyze the case of a rebit ring, we will find that the optimal states, 
which are the states most relevant to our problem, are in fact rotationally 
invariant, so that the above interpretation applies. 
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5 Analysis of a rebit ring 



So far we have not made any assumptions about the structure of the lattice. 
In order to obtain a concrete result, we now specialize to the simplest possible 
lattice for which our general question can be addressed, namely, a closed 
one-dimensional ring. Let the ring consist of n rebits with n > 2, labeled by 
j = 0, 1, . . ., n — 1; the labeling is mod n, so that j = n is the same as j = 0. 
Let Cj be the matrix element c of Eq. (|l^) when the two particles in question 
are particles j and j + 1. The quantity K = ^Y,j \ c j\ wm be the measure 
of average nearest-neighbor correlation that we will be trying to maximize; 
note that < K < 1. The product ]X,( c j/l c j|)) which we will call e ld , is 
the net phase factor associated with transport around the whole ring, and 9 
(defined only mod 2tt) is the net rotation angle. Our question is this: what 
is the maximum possible value of K for a fixed value of 91 Let us call this 
maximum value K ms ^(9). If K max (9) decreases as e td gets farther from unity, 
then we can say that there is a correlation cost associated with non-trivial 
parallel transport. 

As we have stated the problem so far, the phases of the different c/s, that 
is, the phases that define the individual parallel transporters, need not be 
the same for all links in the ring. However, for our purpose there is no loss 
of generality in assuming that these phases are all equal. This is because if 
they were not equal, we could always apply local rotations to the individual 
rebits (analogous to a gauge transformation) so as to make them equal. Local 
rotations can change neither the magnitude of any Cj nor the overall phase 
factor e l6 . So the restriction to equal phases does not eliminate any states 
that might change the answer to our question. That is, states with maximal 
K for any given 9 are still represented in the restricted set. 

With this restriction, we can simplify our problem by expressing the av- 
erage correlation K and the overall phase 9 in terms of creation and annihi- 
lation operators. Let aj be an operator on particle j defined by Oj|t) HI) 
and <ij\l) = 0. Then for the link between the jth and the (j + l)st rebits, we 
can write Cj as 

cj = Tripa^aj), (23) 
p being the density matrix of the ring. Now if we define the operator V to be 

r = ;EsU ( 24 ) 

3=0 
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and write (r) = Tr pT, it follows (assuming that each Cj has the same phase) 
that 

K = 2\(T)\, (25) 
and that the overall phase 9 (mod 2ir) is simply 



9 = narg((r)). 



(26) 



T is not a Hermitian operator — so its eigenvalues may be complex and indeed 
must be complex if 9 is to take any non-trivial value — but V does commute 
with its adjoint, which implies that eigenvectors corresponding to distinct 
eigenvalues are orthogonal. Note also that V commutes with the total z 
component of spin; so the eigenstates of V can be taken to be states with a 
definite number u of up spins. Our immediate goal is to find the eigenvalues 
of r, from which we will be able to determine the set of possible values of 
(T), which in turn will give us K milx (9). 

The a operators are not quite fermionic creation and annihilation opera- 
tors, because the operators associated with different sites commute with each 



other rather than anticommute. However, we can use a standard trick |2T 
to define genuinely fermionic operators bf 



cxp 



i7T E a k a k 
k=0 



(27) 



The operators bj satisfy the usual fermionic anticommutation relations: 

{b J ,b k } = {b],b{}=0; {&* = (28) 

We now express T in terms of the b's. The expression depends on whether 
u, the number of up spins in the ring, is even or odd; that is, the expression 
is different in different subspaces. For odd values of u, T looks the same in 
terms of the fe's as it does in terms of the a's: 



n-l 



(29) 



j=0 

For even values of u, there is a sign change in the last term: 

n-2 



E h Uh - b\b 

j=0 



n-l 



(30) 
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In either case, we can diagonalize T and find its exact single-fermion eigen- 
values. From these we can obtain the eigenvalues of V for an arbitrary value 
of u by summing u of the single-fermion eigenvalues. 

For odd u, the single-fermion eigenstates of T are \uj m ) = rfJjO), m = 
0, ... ,n — 1, where |0) is the vacuum state, that is, the state with all spins 
down, and the creation operator d} m is given by 

dl = ^£e- 2 ™ ij,n b). (31) 

v n j=0 

The corresponding eigenvalues of T are 

g m = {l/n)e 2m ™' n . (32) 
For even u, the creation operators for the single-fermion eigenstates are 

i n— 1 

dl = —J2 e- {2m+1 ^'% m = 0, . . . , n - 1, (33) 

and the corresponding eigenvalues are 

g m = (l/n)e^ 2m+1 ^ n . (34) 

Regardless of the value of u, the single-fermion eigenvalues of T are complex 
numbers of length 1/n, with phases uniformly spaced around the complex 
plane. 

The eigenvalues of T corresponding to a system of u fermions (that is, u 
up spins) are all the possible sums of u of the g m 's. That is, we can write 
each such eigenvalue as 

G(M) = 9m, (35) 

meM 

where M is a set of exactly u integers chosen from the set {0, . . . , n — 1}. (We 
cannot use the same value of m twice in this sum because no two identical 
fermions can be in the same state, and all the single-fermion eigenvalues are 
non-degenerate.) 

As we will see, the most important eigenvalues G(M) for our purpose will 
be the ones with the greatest magnitude. For even values of n, these are the 
ones for which u = n/2 and the set M consists of a string of consecutive 



15 



integers (mod n), so that the corresponding values g m constitute a "fan" of 
complex numbers spread out over half of the complex plane. Performing the 
sum in Eq. (|3^), we find that these extreme eigenvalues are 

g2r7ri/n 

G r = — — , r = 0, . . . , n — 1, (even n) (36) 

n sm(7r/n) 

an equation that holds for all even values of n (even though u = n/2 might 
be even or odd). For odd values of n, the eigenvalues with largest magnitude 
are obtained by setting u equal to either (n + l)/2 or {n — l)/2, and again 
letting M consist of a string of consecutive integers mod n. For either of 
these choices of u, the eigenvalues thereby obtained are 

2riri/n 

G r = __cos(7r/(2n)), r = 0,...,n-l. (odd n) (37) 

nsin(7r/n) 

For either even or odd n, each eigenstate corresponding to one of the 
eigenvalues G r can be written as 

where mo = r — \u/2\ mod n. (Here and below, addition in the subscript 
of d is always mod n.) This value of mo places exp(2riri/n) in the center of 
the fan of complex eigenvalues g mo , . . . ,g mo+u -i. For later convenience, we 
define the following density matrices based on these eigenstates. For even n, 

p r = |f2 r)n /2) (^r,n/2 1) (39) 

and for odd n. 



|^r,(n+l)/2)(^r,(n+l)/2| + |^r,(n-l)/2/ (^r,(n-l)/2 | 



(40) 



As we will see shortly, both of these density matrices satisfy the rebit condi- 
tion, even though |fi r ,(n+i)/2) and |fi r ,(n-i)/2) do not. Note also that for both 
even and odd n, Tr p r T = G r . 

One can check that p r is translationally invariant and thus is in accord 
with the assumption we made earlier, that the matrix element Cj has the 
same phase for each link of the ring. But we also want to check that p r is a 
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legitimate rebit state, i.e., that it satisfies Eq. ([□]). Let us do this first for 
even values of n, in which case we are dealing with a pure state \Q r ,n/2)- 
We begin by noting that 

= (-l)Vf 1 , (41) 
which can be seen directly from the definition of bj. It follows that 

(42) 

where mi, . . . , m u are any distinct values chosen from the set {0, . . . , n — 1}. 
(The addition of n/2 in the subscripts on the right-hand side comes from 
( — iy in the preceding equation.) In general, a state of the form • • • d) |0) 
will not satisfy the rebit condition, even if u = n/2. However, for the spe- 
cial case in which mi, . . . , m u are n/2 consecutive integers, as they are in the 
definition of |fi r n / 2 ), the subscripts on the right-hand side of Eq. ( f42| ) are pre- 
cisely those elements of {0, . . . , n— 1} that are not included in {mi, . . . , m u }. 
Therefore, when those annihilation operators are applied to the all-up-spin 
state <j® ti |0), which within a phase factor is the same as 44"--4-i|0), the 
resulting state, again up to an overall phase factor, is <i] Tll • • • dj^jO). We have 
thus shown that 

vTK,n/2) = e^ n/2 ) (43) 

for some phase /3, so that p r satisfies the rebit condition (|10D for even values 
of n. 

Turning now to the case of odd n, one can use an argument like the one 
in the preceding paragraph to show that 

cr x n |^r,(n-l)/2)(^r,(n-l)/2|0"f n = I ^r,(n+l)/2) (^r,(n+l)/2 1 > ( 4 4) 

and vice versa, so that 

°x PrOx = Pv I 45 ) 

So p r satisfies the rebit condition for odd values of n as well. 

We are now in position to find the set — call it Q — of possible values of 
(T), from which we will be able to determine K max (8). The complex numbers 
G r given by Eq. (|36|) or Eq. (|37|), being values of (r) corresponding to the 
legitimate rebit states p r , are elements of Q. By taking mixtures of these 
states, we can obtain other possible values of (r). Let 

P = J2lrPr: (46) 
r 
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a 



9n d^ ■ ■ ■ d^ 

mi m 



|0> 



"mi 



hn/2 ' ' ' d mu 



K»/2*ri°> 



where the g's are non-negative numbers summing to 1. For this state we have 



(D= Tr P r = ]Tg r G r . (47) 

r 

The complex numbers G r are the vertices of a regular n-gon in the complex 
plane, and Eq. fl4"T|) shows that this n-gon and its interior are contained in Q. 

In fact it is easy to see that Q contains no other points. Any complex 
number (r) in Q must be a weighted average of eigenvalues of F: 

(T) = J2qmG(M). (48) 

M 

But one can show that each eigenvalue G(M), regardless of the value of u, 
lies on or inside the n-gon defined by the special eigenvalues discussed above. 
Therefore it is impossible for the average to get outside this region. 

For the special case n — 2, the interior of the "n-gon" is simply a segment 
of the real axis, running from —1/2 to +1/2. Thus the only possible phases 
of (r) are zero and n, so that according to Eq. (|26|) the only possible value 
of exp(i#) is 1. (This is because there is no real loop to traverse; to return to 
the starting place, one has to retrace one's steps.) For all other values of n, 
all values of 9 from to 2ir are possible. To see this, it is enough to consider 
a single side of the n-gon; let us take the side consisting of the line segment 
joining the point Go, on the positive real axis, with the point G\ = e 2m ^ n Go. 
As we travel along this segment, the phase of (r) varies from to 2n/n, so 
that 9 varies from to 2tt. The range of values of | (r) | as a function of 9 is 
the same for each of the other sides of the n-gon. 

It thus becomes a simple geometric problem to find K max _(9). For n > 2, 
consider the line segment just described, connecting G to G\, and note that 
for any 9 in the range < 9 < 2tt, K max (9) is twice the magnitude of the 
unique point along this segment whose phase is 9/n. Doing the geometry, 
and using the values of G r given in Eqs. (|36D and ([37]), one finds that for 
even n, 

K (9) = 2c ° s(7r/n) (49) 

n sm(7r/n) cos[(7r — 9) jn\ 



and for odd n, 



2cos(7r/n) cos(7r/(2n)) 

A max (fl) = — , f — -. (50) 

n sm(7r/n) cos[(7T — 9)/n\ 
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It is clear both from the geometric picture and from Eqs. (f49|) and ( |50| ) that 
ii'max is largest at 9 = and 9 = 2n and smallest at 9 = it. Indeed, the value 
of K max becomes smaller the more e td differs from unity. In this sense there 
is a correlation cost of non-trivial parallel transport around the ring. 

It is not hard to interpret the states p r physically. The state p , which 
entails no twisting as one goes around the ring, is the ground state, or in 
the case of odd n an equal mixture of the two degenerate ground states, of 



the ferromagnetic XY model [21] on a one- dimensional ring, whose Hamil- 



tonian is H = — J2j( a j a j+i + oJ+iOj) = —n(V^ + T). The state p r can be 
obtained from po by rotating each rebit, the rotation angle at site j being 
2%rj/n. These rotations do not change the strength of the nearest-neighbor 
correlations, but for each site j they change the phase of the matrix element 
Cj from zero to 2irr/n. Still, this does not change the overall phase factor e %e 
associated with the whole ring. When one creates a mixture of two of these 
differently rotated states, e.g., p and pi, the resulting matrix element Cj is 
an average of two complex numbers with different phases. It is this averaging 
process that allows the possibility of a non-trivial net phase change around 
the ring. 

In the case of even n, where p r represents the pure state \ft r ,n/2), one 
could achieve the same averaging effect by creating coherent superpositions 
of eigenstates of T rather than incoherent mixtures. I have chosen to use mix- 
tures because superpositions of these eigenstates are not necessarily trans- 
lationally invariant. However, just to demonstrate that it is possible for a 
translationally invariant pure state to have a non-zero rotation angle 9 asso- 
ciated with parallel transport around the ring, I offer the following example 
for n = 6: 

i*> = -^=[e^ 2 (imm> + •■■) + e-^ /2 amm> + •••)]■ (51) 



Here each ellipsis stands for all possible translations of the given state. 
(Though we are not paying particular attention to the overall phase fac- 
tors of pure states, I have chosen the overall phase in Eq. ( |5T| ) to satisfy 
Eq- (§)•) One finds that the matrix element Cj for each link in this ring is 
(1/6) exp(z0), so that the overall rotation angle is 9 = 6<p and the average 
correlation is K — 1/3. Thus any value of 9 can be realized with a trans- 
lationally invariant pure state. But the value K — 1/3 is not optimal. To 
obtain the optimal value K mSuX (9) in a translationally invariant state, one 
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must typically use mixed states rather than pure states. 

We finish this section by giving asymptotic expressions for K mSjX (9) as the 
number of rebits in the ring gets very large. For even n, Eq. ( fj9|) to order 
1/n 2 becomes (for < 9 < 27r) 



m = - 

7T 



l + -f- 



n- 



while for odd n we have 

^max(^) 



1 + 



6 



7T 



7l9 + — 



n 2 V24 



-tt9 



°1 
2 



(52) 



(53) 



Thus the correlation cost of non-trivial parallel transport becomes smaller as 
the size of the ring increases. 



6 Other lattices 

Let us now think about how the above problem might be generalized to a 
finite or infinite lattice of higher dimension (still using rebits as our basic 
objects). We can state the problem as follows. As in Section [|, let p^ k ^ 
be the density matrix of the pair of rebits at the neighboring sites j and k, 
and let Cj^ be the coefficient of |TI)(1T| i n this density matrix. (This Cj t k is 
analogous to the Cj of the preceding section.) Now suppose that the phases 
of all the c's, for all the links {j, k), are specified. We will call this complete 
specification A, since it is analogous to a field configuration in the U(l) gauge 
theory. Given this specification, we have two questions: (i) Is it possible to 
find a lattice state for which the numbers all non-zero and have the 

chosen phases? (One can always find a state in which all the Cj^'s are zero, 
making the state consistent with any phases, but such a state is not very 
interesting.) (ii) What is the maximum possible value of 

^=yEM, (54) 

consistent with the specification A! Here L is the number of links in the 
lattice. In the case of an infinite lattice, K can be defined as a limit over a 
sequence of finite lattices. Let us call the maximum value i^ max (^4). 
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For definiteness let us consider a specific lattice, namely, an infinite square 
lattice in two dimensions. Let us consider first the configuration Ao in which 
the phases of all the c^'s are zero. In this case it is again helpful to invoke 
the Hamiltonian of the ferromagnetic XY model: 

H = — ^2 ( a l a j + aja&y (55) 

Here the sum is over all links in the square lattice. The optimal value of K is 
the infinite-lattice limit of (—Eq/L), Eq being the minimum eigenvalue of this 
Hamiltonian. In the thermodynamic limit, the ground state of the XY model 
on a square lattice breaks the SO (2) symmetry of the problem and picks out 
a preferred direction of magnetization in the x-y plane |22|, p3| , which can 
be characterized by a single angle a. But if we choose to do so (in order 
to simplify the interpretation of parallel transport), we can easily generate 
a rotationally invariant state with the same energy — or in our context, with 
the same degree of correlation — simply by averaging the ground state density 
matrix p(a) over all angles: p = (1/ 27r) / p(a)da. The ground-state energy 
has been evaluated numerically |24|, |25], [2~6fl , and one finds that K max (Ao) = 
0.549. Notice that this value is smaller than the corresponding value for a 
rebit ring in the limit n — > oo (see Eqs. ( |5^ ) and fl53|) or Ref. [[?!]]), which is 
2/tt = 0.637. 

Given a different specification of the phases, it is not immediately obvious 
whether there exists a state that has all non-zero c^'s — let us call such a 
state "fully connected" — and that is also consistent with the given phases. 
Consider, for example, the configuration A in which all the phases are zero 
except at a specific link (j, k), where the phase is required to be 0. Can one 
find a fully connected state of the lattice consistent with these phases? The 
following method will work, though it is not likely to be optimal. Start with 
the state po defined in the preceding paragraph, in which all the phases are 
zero. Construct the following sum: 

p(j, k, 0) = \ \V jPo V} + vipoV k + [V] VjpoM ® Vi)] , (56) 



where Vj is the matrix 
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applied to particle j and Vk is the same matrix applied to particle k. Here 
£ will be a function of <fi to be determined later. Because the rotations in 
Eq. affect only particles j and fc, all links not involving either of these 
particles will continue to have zero phase. Moreover, any link involving only 
one of the two special sites will likewise have its phase unchanged. Consider, 
for example, the link (j, I) where / 7^ k. The value of Cji is 

c jtl = i(cbe-« + co + c e*) = ( 1 + 2 ™® )co, (58) 

so that the phase has not been affected. On the other hand, the value of c 
associated with the link (j, k) is 

/2e~* + e 2 *\ 
c 3,k=y 3 Jc , (59) 

which has a phase that can be made equal to (f) by a proper choice of the 
value of £. 

An even simpler strategy, which is surely not optimal, shows that any 
phase configuration A can be realized in a fully connected quantum state of 
the lattice. Let us imagine the two dimensions of the lattice to be horizontal 
and vertical. Start with a state in which each vertical column of lattice sites 
is in the ground state of the ferromagnetic XY model for an infinite chain. 
Call this state \4>y)- Now rotate each of the rebits in each of these chains 
so as to achieve the desired phases for the vertical links. This can be done 
with no loss of correlation, because we are simply performing local rotations. 
Similarly, consider the state \4>h) in which each horizontal row is in the XY 
ground state, and rotate the rebits so as to achieve the desired phases for 
the horizontal links. Let \(p' v ) and \<p' H ) be the states resulting from these 
rotations. Then the mixed state 

P=\{WvWv\ + \<t>'H)(<t>' H \) (60) 

completely matches the phase configuration A. We can even compute the 
value of K for the state p: it is equal to half of K max for the infinite chain, 
independent of the configuration A. That is, K = 1/tt — 0.318. For the 
special configuration considered above, in which only one link has non-zero 
phase, this value is smaller than what one can achieve with the specialized 
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method of Eq. (pTp. Nevertheless, the method we have just described does 
answer our first question: all configurations A can be achieved without mak- 
ing any Cj^ vanish. Notice also that this construction gives us a lower bound 
on K max (A) for all configurations A: K max (A) > 1/tt. 

Actually finding K max (A), even for simple configurations A, is probably 
a very hard problem. If one uses strategies similar in spirit to the one given 
in Eq. (|57D , then it would seem that the value of K mas _(A) must decrease in 
order to achieve non-trivial parallel transport around a loop. (One can see 
the decrease in the values of \cj t i\ and |cj ( fc| in Eqs. ( |58|) and fl59f).) But it is 



conceivable that a completely different strategy could do much better; so we 
must leave this basic question unanswered. 



7 Lattices of qubits 

In our definition of parallel transport for rebits, developed in Section |], we 
implicitly made use of the fact that for rebits, there exists a two-particle 
state p( + ^ such that if particles j and k are in this state, and if particle j 
is measured and found to be in the state \tp), then particle k will always be 
brought to the same state Mathematically, 

^ = ^|P (+) I^> = I^>(# (61) 

The state p<+> is in fact |*(+))(#«|, with |*(+)) = (|t|)+ \U})/V2. So for 
this state, the parallel transporter is the identity and the degree of correlation 
is 1. In effect, our definition of parallel transport compares other two-particle 
states to this special state; when the parallel transporter U is not the identity, 
it is because one of the particles has been rotated (and possibly distorted in 
other ways as well) compared to the standard state p( + \ 

For qubits, there is no two-particle state with this property. The most 
closely analogous state is the singlet state l^ - )) = (|T!) — llt))/\/2- It 
has the property that if a measurement on one of the particles brings it 
to the state \ip), the other particle will be brought to the orthogonal state 
\ip) = o- y \ip*), where the complex conjugation is in the standard basis. In 
defining parallel transport for qubits, we will take the singlet state as our 
standard state, for which the parallel transporter is defined to be the identity. 
All other states will then be compared with this one. A general parallel 
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transporter will be a rotation of the Bloch sphere, that is, an element of 
SO (3); as in the case of rebits, our definition will not allow us to distinguish 
overall phases. We will usually represent such a rotation as a 2 x 2 unitary 
matrix, keeping in mind that the overall phase is irrelevant. 
We start again with Eq. (|14"D : 

«,(*) = I(^0')| p 0- fc )|^0')), (62) 

and again let M. be the map that takes \ip&) to u^ k \ But now we define our 
generalized fidelity as 

F{R) = jp# ' (63) 

where dip refers to the uniform measure over the surface of the Bloch sphere. 
Note that F compares A4(ip) with a rotated version of rather than a 
rotated version of If there is a unique rotation R = U (up to an overall 
phase) that maximizes F(R), then we will take this U to be the parallel 
transporter for the link {j, k) . Thus, if p^ k > happens to be the singlet state, 
we have U(k,j) = I. 

Carrying out the integrals in Eq. (|63"D, we find that 



F(R) = \ + ^ <*<->| (I ® I$)pM(I ® R)\¥^). (64) 

Now, it is a fact that any maximally entangled state of two qubits can be 
written as (I ® -R)!^/^ - - 1 ). So maximizing F(R) over all rotations R is the 
same as finding the maximally entangled state that has the greatest overlap 
with jjw'*'. The quantity 

/ = max(^ ( - ) |(J ® R?)pV k Xl ® ( 65 ) 

has been called the "fully entangled fraction" of p (jfc) |TJ. It ranges from 1/4 
(for the completely mixed state) to 1 (for a maximally entangled state). f\ We 
will take (4/ — l)/3, which ranges from to 1, as our measure of the degree 
of correlation between particles j and k, and define K q (q for "qubit") to 

2 Even though the word "entangled" appears in the description of /, it is not a proper 
measure of entanglement; a separable state can have a value of / greater than 1/4. 
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be the average of this quantity over all the links of the lattice. As has been 
mentioned, the rotation R that achieves this maximum value, if it is unique, 
will be our parallel transporter U(k,j)^\ We would like to maximize K q for 
a fixed set of parallel transporters. 

Our notion of parallel transport for qubits is particularly simple if p^^ 
happens to be a "twisted Werner state," that is, a state of the form [^7[ 



p(/®l/)|* ( - ) )(^ ( - ) |(/®l/ t ) + (l-p)(J/4). (66) 

Here < p < 1, |\J/( - )) is the singlet state, / is the 4x4 identity matrix, 
and V is a unitary transformation acting on particle k. If p^ k > is of this 
form, then the parallel transporter U(k,j) works out, not surprisingly, to be 
the transformation V. Moreover, the weight p appearing in Eq. (|66[) is none 
other than our measure of correlation (4/ — l)/3. Thus both the parallel 
transporter and the degree of correlation are particularly easy to interpret 
in this case. One can show that the six-rebit state of Eq. (fy]), reinterpreted 
as the state of a six-qubit ring, has the property that each pair of nearest 
neighbors is of the twisted Werner form. 

A more interesting example of a qubit ring exhibiting non-trivial parallel 
transport — but whose pairs are not necessarily of the twisted Werner form — 
is given by the following state of six qubits. 

IV) = a(|TTUU> -■•■)+ P[e~*(\ UTUT) -••■) + e i? (|mm> -•■■)] (67) 

+7(imm>- iuuu». 

Here each ellipsis indicates all the translations of the given state, but with 
alternating signs. For example, the coefficient a multiplies 

ITTTUD- ITUUTH ITUITT)- IIUTTTH IUTTTI)- IITTTU)- (68) 

The coefficients a, (3 and 7 are real and positive; their values will be specified 
shortly. Note that each pair of nearest neighbors in this state has the same 
density matrix, so that each has the same degree of correlation and the 



3 In the case of qubits, in contrast to that of rebits, it is possible for R not to be unique 
even when the correlation (4/ — l)/3 is not zero. For example, if p^ fc ) = (|U)(U| + 
|IT)(it|)/2) then any rotation around the z axis maximizes F(R), and yet the correlation 
(4/— l)/3 has the value 1/3. However, as in the case of rebits, this sort of ambiguity does 
not cause any difficulties for the problem we are considering. 



25 



same parallel transporter. Because the state \ip) is an eigenstate of S z with 
eigenvalue zero, the density matrix of each pair is of the form (|T9"D with all 
the Xj's equal to zero. (It is convenient to introduce a negative sign in the 
off-diagonal elements since our standard state is now the singlet.) 
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(69) 



Carrying out the trace of over the other particles, one finds that 

a = 2a 2 + 2(3 2 , b = a 2 + 4[3 2 + 7 2 , and c = 2(3[(a + 7)e*« + (3e~ 2 % The 



correlation K q is (4/ — l)/3, where 

/ = b + \c\ = a 2 + 4(3 2 + 7 2 + 2(3\{a + 7 )e* + (3e~ 2 ^ 
and the parallel transporter is 

U = 




(70) 



(71) 



where is the phase of c, that is, the phase of (a + j)e 1 ^ + (3e~ 2% ^. The fact 
that K q depends on the matrix element b, while the analogous quantity K 
for rebits depended only on c, ultimately comes from the fact that in Eq. (|63|) 
we average over the entire surface of the Bloch sphere and not just over the 
equator. 

In the spirit of our main problem, we would like to choose a, (3 and 7 
so as to maximize K q . Let us do this extremization for the special case 
£ = 0, for which the parallel transporter given in Eq. (|71~1) is the identity. 
In this case one finds that the optimal values are a = (130 + 34-\/l3)~ 1 / 2 , 
(3 = (1/2) (3 + \^L3)a, and 7 = (4 + ^/l3)a, which satisfy the normalization 
condition 6a 2 + 12/3 2 + 2 7 2 = lfj In what follows we will assume that a, [3 
and 7 have these values. One finds then that for £ = 0, each pair of nearest 
neighbors is in a Werner state, and the correlation K q is (2 + V^13)/9 = 0.623. 
Let us call this value K^. 



4 By no accident, the state \ifj) with these values of the coefficients is the ground state 
of the antiferromagnetic Heisenberg model for a ring of six qubits. 
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How do K q and U for the state change as £ departs from zero? Let 
us first look at U. To lowest order in £, the angle 0, which is the rotation 
angle associated with parallel transport across a link, is 

= arg[(a + 7 )(1 + z£) + /3(1 - 2i0] = f " + 7 ~ 2 f V- ( 72 ) 



+ 7 + /? / 



Note that a +7 — 2/? > 0, so that this linear contribution to does not vanish. 
Meanwhile, the correlation K q diminishes by an amount proportional to the 
square of £: 



Kq = K (0) f'12/? 2 (a + 7) 



)e 2 - (73) 



a + 7 + /3 

Letting = 60 be the net rotation associated with parallel transport around 
the whole ring, we can see how K q depends on 9 to lowest order: 

Notice that the reduction in correlation is of second order in 6, whereas in 
the case of rebits it is of first order as seen in Eqs. ([52] ) and (|53|) . Of course 
we have not done the thorough optimization for qubits that we have done for 
rebit rings, but this example indicates that there is a significant difference 
between the two cases. 



8 Discussion 

We have shown, first of all, that there is a sense in which certain quantum 
states exhibit non-trivial parallel transport around a closed loop, which is 
to say that the nearest-neighbor correlations exhibit a net twist as one goes 
around the loop. One might regard this result as somewhat surprising, since 
there is, after all, only a single quantum state for the whole loop, and one 
might think that the local twists would therefore have to cancel each other 
out. We have also shown that in the case of the rebit ring, there is a loss 
of nearest-neighbor correlation associated with non-trivial parallel transport 
around the ring. For other lattices or for lattices of qubits, we do not know 
whether there is such a correlation cost, though it is certainly plausible that 
there would be. 
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As we have seen, there is a close relationship between our rebit problem 
and the XY model, a model that has been very well studied. Studies of spin 
stiffness in this model (see, for example, Ref. |26|| ) have a certain similarity 
with the problem we have been considering in that in both cases one enforces 
a twist between neighboring spins. What distinguishes the present work is the 
fact that we have not actually specified any Hamiltonian. Though it has been 
helpful for us to use an operator similar to the Hamiltonian of the XY model, 
we are really working with what might be called the kinematics of quantum 
mechanics. We are asking what correlation properties of quantum states 
follow from certain other correlation properties — specifically, we are asking 
what one can say about the strength of correlations given some information 
about the twist in the correlations — and this question is independent of any 
considerations of energy. 

I introduced the subject by relating it to lattice gauge theory. To what 
extent, then, have we found an analogy between the state of a quantum 
lattice and the configuration of a lattice gauge field? In a qualitative sense, 
the reduction in correlation that we have observed in a rebit ring can be 
compared to the reduction in probability that one finds in a lattice gauge 
theory. But for this rebit case, the analogy must be regarded as quite rough, 
because there is a significant lack of congruence in the details. In gauge 
theory, the initial decrease in the probability is of second order in the net 
rotation angle 9 associated with transport around a plaquette. In the U(l) 
theory, for example, the function exp(— 5"), with S given by Eq. (H), decreases 
in proportion to 9 2 for small values of 9. In contrast, in the rebit ring, the 
average correlation K, as given by Eqs. ([49D and (|50|) , decreases in proportion 
to 9 itself. The second-order dependence is in fact important in lattice gauge 
theory for taking the continuum limit. Moreover, our first-order dependence 
makes K ma ^(9) a non- analytic function, since near 9 = it takes the form 
a — b\9\. So this difference is not trivial. 

On the other hand, we have just seen that the qubit correlation as we 
have defined it does seem to diminish quadratically in 9, at least for a ring of 
six qubits. It is interesting to ask whether in the case of a two-dimensional 
or higher- dimensional lattice, the dependence of K m3iX (A) is of second order 
in the rotation angles. In the one relevant example we have considered for a 
rebit lattice, namely, the strategy given in Eq. (|57|), the value of K decreases 
as 2//3 , which is an even sharper dependence than in the rebit ring. But 
we have not explored at all fully the range of possible states that one might 



28 



consider for these higher-dimensional lattices. 

It is worth commenting on the fact that in the case of qubits, our approach 
makes the parallel transporters elements of SO (3), whereas one might have 
expected SU(2). The nature of our definition does not allow us to pick out a 
relative phase in the relation between neighboring qubits. For example, we 
cannot distinguish between the identity operation and a rotation by 2n, even 
though a pure qubit state experiencing the latter rotation picks up a phase 
factor of -1. It is conceivable that by taking into account the density matrix 
of an entire loop, in addition to the density matrices of the neighboring pairs, 
one might be able to make sense of this distinction as it applies to the net 
rotation associated with the loop as a whole. 

There are other ways in which one might modify the problem we have 
been considering. One could use a different measure of correlation or entan- 
glement. Moreover, even if one continues to use the quantity 2|c| for rebits 
and the quantity (4/ — l)/3 for qubits as the measure of nearest-neighbor 
correlation, one could combine the correlations from all the links in a dif- 
ferent way. For example, in the case of a ring it would make some sense to 
consider the product of the individual correlations rather that the average; 
this measure has the pleasing feature that it vanishes if any of the links in 
the ring is broken. 

Again, our main conclusion is this: in the one example we have worked 
out in detail, if there is a non-trivial twist in the nearest- neighbor corre- 
lations, that is, a twisting that cannot be undone by local rotations, then 
there is a corresponding reduction in the maximum possible magnitude of 
these correlations. That is, in this one example at least, twisted correlations 
are weaker correlations. This conclusion follows from the structure of the 
quantum state space and is true irrespective of the system's Hamiltonian. 
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